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Abstract 

In the framework of the grand-canonical ensemble of statistical me- 
chanics, we give an exact diagrammatic representation of the density 
profiles in a classical multicomponent plasma near a dielectric wall. 
By a reorganization of Mayer diagrams for the fugacity expansions of 
the densities, we exhibit how the long-range of both the self-energy 
and pair interaction are exponentially screened at large distances from 
the wall. However, the self-energy due to Coulomb interaction with 
images still diverges in the vicinity of the dielectric wall and the vari- 
ation of the density is drastically different at short or large distances 
from the wall. This variation is involved in the inhomogeneous Debye- 
Hiickel equation obeyed by the screened pair potential. Then the main 
difficulty lies in the determination of the latter potential at every dis- 
tance. We solve this problem by devising a systematic expansion with 
respect to the ratio of the fundamental length scales involved in the 
two coulombic effects at stake. (The application of this method to a 
plasma confined between two ideally conducting plates and to a quan- 
tum plasma will be presented elsewhere). As a result we derive the 
exact analytical perturbative expressions for the density profiles up to 
first order in the coupling between charges. The mean-field approach 
displayed in Paper I is then justified. 

KEYWORDS : Coulomb interactions, dielectric wall, grand-ca- 
nonical ensemble, systematic resummations, inhomogeneous Debye- 
Hiickel equation, screened potential with two characteristic length sca- 
les. 
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1 Introduction 



The present paper is devoted to the systematic derivation of the density 
profiles which are discussed in Paper I. The system is a classical multicom- 
ponent plasma, made of at least two species of moving charges with opposite 
signs, near a plane wall macroscopically characterized by a dielectric con- 
stant. The exact analytic expressions are obtained in the limit of weak 
Coulomb coupling inside the fluid. 

Our calculations performed in the framework of statistical mechanics 
cast a new light on the fundamental phenomenon in such systems : Coulomb 
screening of surface polarization charges inside the fluid. Indeed, from elec- 
trostatics it is well known that the bulk thermodynamic properties in a 
fluid of charges are independent from the charge state of boundary walls. 
Nevertheless the microscopic long-ranged Coulomb pair interaction f (r; r') 
between two unit charges located respectively at r and r' takes into ac- 
count the electrostatic response of the wall : f (r; r') is a solution of Poisson 
equation 

Ar'v{r; r') = -47r5(r - r') (1.1) 

and its expression is ruled by the electrostatic boundary conditions. For 
instance, when the material of the wall is characterized by a relative dielec- 
tric constant (with respect to the dielectric constant in the vacuum), the 
solution of ( p!T| ) reads 

vUr; r') = - Aei^— (1.2) 

|r — r'l |r — r | 

where Aei = (£„ — l)/(ew + 1) and r'* is the image of r' with respect to the 
plane interface |^]. On the contrary, far away from any boundary or when 
the wall has no electrostatic response (ew = 1) the potential takes its "bulk" 
value, 

v,{r;r') = -^ (1.3) 
|r — r'l 

Short-distance cut-offs must be introduced in order to prevent the collapse 
between charges with opposite signs or between every charge and its image 
when Ew > 1- 

At the inverse temperature f3 = l/k^T, where kg is Boltzmann constant, 
the four length scales in the problem are : the closest approach distance /3e^ 
determined by Coulomb interaction between charges of typical magnitude e 
and with a mean kinetic energy of order 1//3; the mean interparticle distance 
a; the range b of the hard-core repulsion from the wall (which is chosen to be 
the same for all species); the hard-core diameter a of charges which prevents 
the electrostatic collapse inside the fluid. However the cut-off a proves not 
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to arise in the densities and correlations at leading order in a low-density 
regime (because the long range of Coulomb interactions is then the most 
important effect). 

In the grand canonical ensemble (Section the macroscopic parameters 
are the volume of the system, the fugacities Za's, where a is a species index 
running from 1 to the number of species n^, and the inverse temperature 
/?. The task of getting exact analytic results for a multicomponent plasma 
where a cut-off is added, and where the screening of the long-range part 
of the interaction is systematically dealt with, was achieved for the density 
expansion of bulk correlations by Meeron |11]. Haga extended Meeron's 



scheme to the low-density expansion of the pressure [g] in the bulk where 
densities are uniform]]. 

In Section ^ we devise a new reorganization of Mayer fugacity-expan- 
sions for every particle density /3„(x; {z-y}, /3) and its uniform bulk limit 
Paii^'y}, (3) in order to exactly deal with the long range of the self-energy 
of a charge interacting with its image inside the wall as well as with the 
long range of Coulomb pair interactions, (x is the distance from the plane 
wall located at x = 0. The notations {z^} means that the fugacities of all 
species are involved.) Our reorganization of diagrams is performed in two 
steps. In the first step we exhibit how the screening of the bare one-body 
interaction with the dielectric wall may be described by resummations of 
"ring" subdiagrams. In the second step, we use an already known method 
in order to exactly handle with the screening of Coulomb pair interactions 
[§]. (We notice that, if we had used a one-step resummation analogous to 
that performed in Ref.[|[| or Q, there would have been resummed diagramms 
with some bare self-energies left - and therefore some spurious 1 /x tails left- 
which would disappear only by an adequate gathering with other diagrams 
to be specified in every particular case.) 

The latter exact systematic resummations produce auxiliary screened 
potentials 4>i and (p2 each of which obeys a second-order differential equation 
studied in Section ^. The technical difficulty lies in the resolution of these 
auxiliary inhomogeneous Debye-Hiickel equations. 

In order to get explicit analytic expressions, we consider a regime of weak 
coupling inside the fluid 

r = ^ < 1. (1.4) 

a 

However, the electrostatic coupling with the dielectric wall is not necessarily 
weak even in the very vicinity of the wall, 

^ - 1 (1-5) 



^ Apart from a minor error in the calculation of the coefficient called S2, his result is 
correct. 
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After resummation of coulombic long-range divergencies, collective screening 
effects operate on a scale equal to the bulk Debye length, k~^^, with 



1 

r2 



oc 



(1.6) 



The weak-coupling condition ( |1.4| ) is equivalent to the condition 



(1.7) 



At the same time, the density profile proves to be an expansion in integer 
powers - possibly multiplied by some powers of logarithms - not of the 
coupling parameter /3e^/a = F but of 



(1.^ 



Thus an e^-expansion provides an expansion in powers of the square root of 
the density. (The cut-off a for pair interactions corresponds to an integrable 
interaction ; so it arises in the corrections to bulk quantities only from the 
relative order proportional to the density, namely from the order .) 

We show that the screened potentials 4>i and (p2 can be expanded in 
powers of the ratio A of the two length scales l3e'^ and which controls the 
variations of the x-dependent inverse Debye length in the equation obeyed 
by each The expansion of each (f)i(r,r') is uniform with respect to the 
distances x and x' from the wall - when the projection of r — r' onto the 
wall is kept fixed - in the sense that the term of order A" is bounded for all 
X > b and x' > 6 by a constant which depends on Kj^b and X/[Kiyb] = /3e^/6. 
In fact A coincides with Ed in the present problem and a scaling analysis of 
Mayer diagrams allows one to obtain the small-So expansion of the density 
profile (see Section |5|). Only one diagram contributes up to first order in 
and paix; {z-y} , 13) can be explicitly rewritten as pa{x; {p^} , (3) . The latter 
general scheme, which provides systematic expansions in A and Ed, can be 
also applied to a classical plasma near a charged dielectric or an ideally 
conducting wall or to a low-degenerated quantum plasma in the vicinity of 
a boundary with any electrostatic response. The corresponding works will 
be displayed elsewhere. 

Eventually, at first order in the coupling parameter e^,, 



Pa{x) = pIO{x - b) exp 



2k£)X ' 



Ax 



X [1 + R{KDx;eD,KDb)] (1.9) 



where R is of order eolneo and e^. R{kox) decays exponentially fast 
at large distances over the scale and tends to a finite value when x 
comes to zero. Subsequently, paix) varies drastically over the bare-coupling 
scale /3e^ very close to the wall, whereas its variation at distances larger 
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than is scaled by the screening length k^^- In other words, very near 
the wall, paix) is governed by the part of the self-energy of each particle 
originating from the electrostatic response of the dielectric wall, whereas the 
mean Coulomb interaction with other charges of the fluid modified by the 
impenetrable wall dominates the approach to the bulk value far away from 
the wall. After calculation of the electrostatic potential drop ^(x) at leading 
order e^, from the charge density profiles given by (1.9), the density profiles 
prove to coincide with the expressions derived from the mean-field approach 
introduced in Paper I, 

p^(x) = pI e{x - b) e-'3'='KdfW [1 _ /3e„$(x)] (1.10) 

where V"^^^{x) is the screened self-energy at order arising from the first 
step in the resummations of Mayer diagrams. 



2 Model 

2.1 Potential energy 

The total potential energy is the sum of the electrostatic energy f/giect and 
the short-range repulsive interaction Us-r- U^iect is equal to the sum of 
Coulomb pair interactions plus the sum of self-energies in the presence of 
the wall, 

^^eiect = ^ ^eQ,,e«^.t;w(ri;rj) + ^ e^^Fseif (a;^) (2.1) 
In ( p.l[) Cq,- is the charge of particle i and 

Kelf(xi) = -Ael^ (2.2) 

where Agi is defined in (|1.2D . (see Ref. 1^] ) . The repulsion energy Uq^ is also 
the sum of hard-core pair interactions plus the sum of one-body potentials 
created by the wall, 

^sk = \Y. ^sR(ri - rj) + ^sR(a;^) (2-3) 

i^j i 

with 

exp[-/3...fe.,)i={° (2.4) 

and 

exp [-pVs^ix)] = I J ![ ^ < ^ (2.5) 

We notice that in the case ew < 1, exp [—Pe^Vs^if^x]) vanishes when x goes 
to zero and Vsr can be omitted, i.e. b can be set to zero, without any collapse 
of the system onto the wall. 
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2.2 Fugacity expansions 



The grand partition function of the system in the semi-infinite space x > 
reads 



+00 

N=0 ' 



N 



JJ dVi Za^ 



.i=l 



-/3(C/clcct + (^SR) 



(2.6) 



p ria 

In (^]^) the integration measure is denoted by / dVi = / dri and 

Jx,>0 ^^^^ 

Vi is the notation for (rj,aj). z^. has the dimension of an inverse cubed 
distance. The one-body interactions can be absorbed in an x-dependent 
fugacity 



Zrv^ \ j Zfy- 6 



-l3[Vsii{xi)+el.V^^l{(x,)] 



According to and (|2 



ZaiO{xi - b)e 



(2.7) 



(2i 



where 9 is the Heaviside function. H can be written as 



+00 

N=0 ' 



N 



.1=1 



_§_ 

e 2 



(2.9) 



The particle density for species a is derived from the grand partition function 
through a functional derivative 



Pa{x;(3,{z.,]) = Za{x) 



5 In H (A {4}) 



(2.10) 



6z*{x) 

By definition of the bulk 

lirn p^{x;(3,{z^}) = pI{P,{z^}) (2.11) 

According to the survey of the existence of thermodynamic limit for Coulomb 
fluids [10|, the bulk densities satisfy the local neutrality 

^ec.pl{(3,{z^})=0 (2.12) 

a 

whatever the values of the In the following we will consider a 

weak coupling regime where the system behaves as an ideal gas at zeroth 
order in T. Besides, the particle density of species a in an ideal gas is just 
proportional to the fugacity Za. As a consequence, we can take advantage 
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of the freedom of choice for the z^'s exhibited in Ref. |10[ | and enforce the 
extra condition 

"]e^z^ = (2.13) 



( 2.13| ) ensures that the ideal gas also satisfies ( 2.12 ) and this provides a 



convenient simplification in the relation between the p^'s and the z-^'s. We 
stress that the irrelevance of the value taken by one among the fugacities, 
which is always valid for bulk quantities, also holds for surface statistical 
properties but only when the wall gets no global influence charge in the 
presence of the Coulomb fluid. This is indeed the case when e„ is finite. 
On the contrary, when is sent to infinity, the wall material is an ideal 
conductor which gains a global charge by infiuence. The latter is exactly 
compensated by the net surface charge in the fiuid and it depends on all 
fugacities z^'s. 

2.3 Representation of the density profile in terms of Mayer 
diagrams 



The fugacity-expansion ( |2.9| ) may be expressed in terms of Mayer diagrams 
where a point is associated with every Vt and two points are joined by at 
most one bond 

J(p.^-p^.) = g-/3(i'SR(ri-rj)+ec.ec^,i;w(ri;rj)) _ ^ (2.14) 

The relation (|2.10| ) and topological arguments lead to various diagram- 
matic representations for the fugacity expansion of Pa{x). In the following 
we will use 



Pa{X) = Za{X) exp 



1 r ^ 



(2.15) 



In ( 2.15| ) the sum runs over all the unlabeled and topologically different con- 



nected diagrams G with one root point V = {a, r) (which is not integrated 
over) and N internal points {N = 1, • • • , oo). Moreover the G's satisfy the 
extra constraint that the root point V is not an articulation point. (V is an 
articulation point of G if G is split into at least two pieces when the point V 
is removed.) [[{ f](^ is the product of the / bonds in the diagram G and Sq 
is the symmetry factor of G, namely the number of permutations between 
some internal points Vn that do not change the product [H/Jg' 

At large distances v^{ri;rj) behaves essentially as l/|rj —rj\ and the / 
bonds ( p. 14 ) are not integrable. The divergences in Mayer diagrams G may 



be removed by splitting the / bonds into several auxiliary bonds /* and 
by performing a reorganization of the diagrams G* made with /* bonds. 
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( G* diagrams have the same topological properties as G diagrams.) The 
decomposition of / chosen for our purpose is 

f = r + l[rf + fT (2.16) 

with 

r%Vi,Vj) = -Pea^ea,v^{r,;rj) (2.17) 

With these definitions /t is at the border of integrability. The systematic 
partial resummations of diagrams G* are described in next section. 



3 Exact topological resummations valid for any 
density 

3.1 Step 1 : Screening of the self-energy induced by the wall 

The long-range one-body interaction e^V^„if (x) involved in the fugacity Za{x) 
turns into a short-range effective interaction when Coulomb "ring" subdia- 
grams are resummed. We introduce the following definitions. A Coulomb 
"ring" subdiagram carried by a point V is either a bond [/^'^J^ /2 or a chain 
of bonds f^^ whose both ends coincide with V and which contains at least 
two intermediate points. The value 1^ of the sum of all rings attached to a 
point V is equal to 



1 



pel [{(pi - v„) (r; r)] 



(3.1) 



where — /3e^(/)i(r; r') is the sum of chains made of an arbitrary number of 
bonds f^^ between r and r' and with intermediate points Vi weighted by 
'ziVi). (1/2 is the symmetry factor of a single ring diagram.) By definition 



,^i(r;r') = /"(P;P') + 



E 

N=l 



N 



\{dViz{Vi) 



.i=l 



nV;Vi)nVr,V2) ■ ■ ■ r{VN;V') (3.2) 



According to (|2.17[) , the definition (|3.2|) of 4>i is equivalent to the integral 
relation 



,{r;v')=v^{r;r')-(3j2^» J dr'%(x>„(r; r")0i(r"; r') 



(3.3) 



The sum of an arbitrary number of rings carried by a point V is equal to 
exp/i.(r). (An analogous calculation appears in section V.B of Ref.[^].) 

Let us consider all diagrams G* (with weight z) that only differ from a 
given diagram G*''^! by one or more Coulomb ring subdiagrams carried by 
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at least one point of G*''^] (see Fig.|^). Before integration over the internal 
points of G*M, the integral corresponding to such a diagram is equal to the 
integrand associated with G*^^^ times the contributions from Coulomb rings. 
As a consequence, if we sum all such diagrams G* we get one diagram G*I^1 
which has the same points and bonds as G*^'^^ except that the weights z are 
replaced by weights 



(-p) = z{V) exp 



--(5e 



.)(r;r) 



(3.4) 



where (/>i(r; r') is defined in (3^). At the end of the first partial resummation, 
the diagrams G* with bonds J^'^, 1/2 [Z^'^]^ and /x and with weights ^q(x) 
have been replaced by diagrams G*'^! with the same kind of bonds but with 
weights Zq^(x) and an extra construction rule 3i (see Fig.|2|). The latter one 
is necessary to avoid double counting. If some subdiagram S(7^m) of G*[^] is 
a ring of bonds /'^'^ with weights z^' carried by the point Vm (with weight 
^al(2;m))) then after integration over intermediate points in the ring, the 
contribution from S(Pm)j which depends only on Xm, is factorized in the 
total contribution from G*'^] and is multiplied by the weight ^aj^(xm) of 
its root point Vra\ rule 3i states that the contribution from S('Pm) is equal 
to its value where all its intermediate points have a weight 'z^ minus the 
corresponding value where all its intermediate points would have a weight 
Zq-. Indeed the latter value is already taken into account in the weight 
'^m^^m) of the point Vm^ whereas a subdiagram of G* made of bonds /'^'^ 
but where at least one intermediate point carries another Coulomb ring does 
not disappear in the ring resummation process leading to G*[^l. (See Figs.| 
and H) . 

According to the definition (2/7) of ^('P) and since V^^iix) = (1/2) {v^ 



-Vb) (r; r), (3^) may be written as 

b)za exp 



[X 



(r;r) 



(3.5) 



As shown below in a weak-coupling expansion, we expect that for any finite 
coupling constant {4>i — fB)(r;r) has no algebraic tail at large distances. 
Thus, contrary to Za{x) the screened fugacity Zo^(x) does not contain the 
long-range part of the self-energy induced by the electrostatic response of the 
wall when ew 7^ 1. In other words, the resummation of rings subdiagrams 
has captured the main effect of Coulomb screening upon the electrostatic 
interaction with the wall. 



3.2 Step 2 : Screening of pair interactions 

In the second step we sum all chains of bonds /'^'^ in diagrams G*[^L The 
systematic resummation process is performed as follows. We define a "Cou- 
lomb" point as the intermediate point of a chain of two bonds /'^'^ and 
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which is not hnked to any other point in the diagram. Such a point has 
a weight Zq^. A so-called "prototype" diagram P is a diagram G*^'^'^ which 
contains no Coulomb point. We sum all diagrams G*[^l which can be built 
from the same diagram P by addition of at least one Coulomb point (with 
weight zt^l) with the associated /'^'^ bonds either to replace a f^'^ bond 
in P by a chain of f'^'^ bonds or to multiply a bond in P (see Fig. 3). In 
the latter case, we use the convention according to which the diagram P 
is the same whether the diagrams G*'^] contain a bond (1/2) [Z'^'^]^ ('^ij'^j) 
or a bond fTi'Pi,'Pj) which are possibly multiplied by a chain of f'^'^ bonds 
between Vi and Vj. Since the P diagram may contain articulation points, 
the resummation produces two effects. 

First, we sum all diagrams G*[^l that differ from P by a ring of bonds f'^'^ 
with at least one intermediate point with weight z'^l. As in the reorganiza- 
tion of diagrams performed in step 1, this process leads to a renormalization 
of the fugacity by an exponential factor analogous to ( |3.4| ) with an extra 
subtraction arising from the construction rule 01 for G*[^l diagrams. The 
resummed fugacity reads 



z^a\x) = exp 



r:r 



v^){r;r)]} (3.6) 



where (t)2 is defined by the integral equation ( |3.3|) with Za{x) replaced by 
Za'(x). According to ( p. 4] ) and (2^), the expression of ^^'(x) is similar to 

TBI) 



b)za exp 



h - WB)(r;r) 



(3.7) 



The second effect of step-2 resummations is a renormalization of the 
/* bonds in the final diagrams P. The summation operates for each pair 
of points in P independently. (See for instance 0] for similar topological 
considerations.) The sum of all possible single chains made of bonds 
and whose intermediate points carry weights z^'^'^ is 



(3.8) 



The sum of the bond (1/2) {V; V) and of all subdia grams made of the 

product of a Coulomb chain (possibly made of only one bond f^^iV'jV')) 
and a Coulomb chain with at least one intermediate point linking V to V' 
is merely 



1 



[F'^'f {V;V') 



(3.9) 



The sum of the bond /t('P;'P')) its product by at least one Coulomb chain 
with at least one intermediate point, and the product of at least three 
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Coulomb chains, one of which at least contains one intermediate point leads 



to the simple resummed bond (see Ref. [11| and Q ) 

... 1 



1 



(3.10) 



After resummations the diagrams P are made of bonds F''^'^, (1/2) [F'^'^]^, 
and with at most one bond between two points (see Fig.^) and the 
following rules. A point - different from the root point - which carries 
only two bonds F''^^ has a weight zU — Za' (3^1 ) (see Fig.|). A point - 
different from the root point - which is linked only to a (1/2) [F'^'^]^ bond 

has also a weight 'z^ — (3^2) (see Fig.^). The rules 3?i and 3^2 arise from 
the definition of the Coulomb points which disappear in the resummation 
process and from rule 3? of the first step of the resummations. All other 

—[21 

points have a weight Za . For instance, the diagrammatic representation of 
/9Q,(r) in terms of P diagrams starts as 



P«(r)=^,^'' 



f (:E)exp|^ y dv'^^^\x')\F--^F^^\ {V;V') + 



dv' 4^\x')-z^^\x') ^[F-(P;P')] +••• (3.11) 



The diagrams represented by the dots contain more than one internal point. 



4 Screened potentials 

4.1 Partial derivative equations 

The screened potential 02) as well as the other auxiliary object (j)i, are 
defined through integral equations (p. 3D. Since is a solution of Poisson 
equation (|1.1|), ( |3.3D is equivalent to a set of local equations. For x' > 0, the 
latter ones read 

Ar(/)i(r; r') - Kl^(x)(/>i(r; r') = -47r(5(r - r') if x > 

(4.1) 

Ar(/)i(r;r') = if x < 

with 

W\x) = 47T/3Y,elza{x) (4.2) 

a 

According to its definition {^^), 4>i is real. Since the operator [Ar — Kr^(x)] 
is self-adjoint, the real function (/>i(r; r') is symmetric under exchange of its 
arguments when r and r' are in the same region, 

</.i(r;r') = <Ai(r';r) (4.3) 
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Moreover, the invariance of the system under translations in directions par- 
ahel to the plane interface implies that 

'Ai(r;r') = 0i(x,x',y - y') (4.4) 

where y is the projection of r onto the plane perpendicular to the x-axis. 
According to its definition ( |3.2| ), 0i obeys the same boundary conditions as 
the electrostatic potential v^: (pi is continuous in all space, 

lim e„^p-{x,x',y -y') = lim ^p-{x, x' ,y - y') (4.5) 

x—fO~ OX x—*0+ OX 

Since ki'^{x) contains a Heaviside distribution at x = b but no Dirac distri- 
bution, the derivative of (j)i{x,x';y — y') is continuous at x = b. Moreover 
4>i{x,x';y — y') tends to zero when x goes to +00 or —00. All previous 
properties are also true for (/>2(r; r') with the only difference 

K^2(^)^4vr/3^e2zW(x) (4.6) 

a 

4.2 Scaling property 

A Fourier transform allows one to change the system of partial derivative 
equations ( [4.1[ ) for each (j)j{r;r') (with j = 1,2) into a system of one- 
dimensional differential equations with respect to x. For x' > 0, 



dx'^ 

dx"^ 
dx^ 



;i + q2)-[/j(x) 



b] {x, x', q) = —4tt5 [x — x') 



i>j {x, x' ,q) = — AttS (x — x') 
lj{x,x',q) = 



ifb < X 

(4.7a) 

if < X < 6 

(4.7b) 

if X < 
(4.7c) 



In ( \l.7\j we have introduced the dimensionless variable with 



/4vr/3^ 



K2 = ATr/J^elzae" 



(4.8) 

(4.9) 
(4.10) 
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and the dimensionless Fourier transform 



(4.11) 



with (j)j{x, x' , Kjq) = J dy exp [iKjq.y] x', y). According to ( p.8D and 
(4.2), for X = Kix 



Ui{x) 



(4.12) 



while, according to (3.5) and (|4.6|), for x = K2X 



TT (~\ '^^^ 2 



^{l/2)l3el{4>l-v^){r;r) 



(4.13) 



When x goes to oo, (/>i tends to the bulk value (pi^ which is the solution of 



[Ar-K?] 0iB(r,r') = -47r5(r-r') 
for all r's. As rederived below, 

0iB(r, r') = Ki0B(«;i|r - r'|) 

with 

e-|f-r'| 



!(f,r') 



r — r 



As a consequence, 

((AiB -^^b) (r;r) = -ki 
and the definitions ( [4.12|) and ( [4. 13]) ensure that 

lim U^{x) = 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



The boundary conditions for 0j(x, x', q) along the x-axis are the same as for 
(/)j(x,x',y — y') (see (^)). Moreover, since K2 = K2{Ki,e) a mere scaling 
analysis of (4.7) shows that 

(pj{r,r';Kj,(3e'^,b,Aci) = Kj(j)j{K.jr, Kjr';e, Kjb, A^) (4.19) 

where (e = (l/2)Ki/3e^) is the generic notation for the c^^'s (sq-'s). 
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4.3 General solutions 



In the following only the value of (/>j(x, x'\ y — y') for x' > b will be involved 
and ( 4.7b| ) becomes an homogeneous equation as well as ( |4.7c| ). The simple 
general solutions of the latter equations that satisfy boundary conditions 
take the form 



(^j(5,5',q; b) 



i?j(x';g;6)(l-Aei)e^* 
Bj{x';q;b) [e'?* + Aeie""*] 



if x < 
if < X < 6 



(4.20) 



where q = |q|. The general solution of (4.7a) is the sum of a particular 
solution (f)"- (x , x' , q) and the general solution hj of the corresponding homo- 
geneous equation 



dx"^ 



x;q) + [Eq-C/j(5;)]/ij(x;q) =0 



(4.21) 



with = — (1 + q^). We look for a solution [h- ) which vanishes (blows 
up) at large positive x, 



lim 



[x;c[) 







(4.22) 



while lima.^_|_oo h^{x; q) = +oo. 

Since the L'j(x)'s vanish when x goes to +oo, then, for any given positive 
number -q > there exists some Xq^^ such that for all x > Xq^rj : Eq — Uj (x) < 
l-Eql+r/. As a consequence (see Ref.|12|) there exists one particular solution 



/ij (x;q) (with a particular multiplicative constant) which tends to zero 
when X goes to +00 at least as fast as exp [— xy^lE'q 



Vx > Xq^^ (x;q)| < exp 



(4.23) 



whereas all other solutions grow to infinity at least as fast as 
exp [xy^|£'q| — rj ]. Let call hj'*{x;q) such a particular solution which di- 
verges when X goes to -|-oo. Henceforth (see for instance ||T3]) a particular 



solution of ( 4.7a ) for x > is just 

4tt 



i*(x,x',q) 



J.q 



hr* (inf(x, x'); q) h^* (sup(x, x'); q) (4.24) 



hr*(3:){dh^*{3:)/dx) - h+*{3:){dhr*{3:)/d3:) is 



72 / J~2 



where the Wronskien Wj^ 

independent from x because d^/dx^ -|- [E^ — {7j(x)] is a self-adjoint operator. 

Finally, a generic solution of ( |4.7aD which vanishes when x goes to -|-oo 
takes the form 

0j(x, x', q; 6j) = (^*(x, x', q) -I- ^(x'; q, 6j)/ij^*(x; q) for 6 < x (4.25) 
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with bj = Kjb. Moreover the symmetry property (4.3) imphes that the 



solution of ( ^.Tal ) with adequate boundary conditions at +00 may be written 
as 

0j(x,i',q;6j) = <A*(i,5',q) +Z(q,6j)/i+*(i;q)/i+*(i';q) (4.26) 
for b < X and b < x'. The coefficient Z(q, b) is determined by the continuity 



relations at 5; = 6j and x = with the solutions (4.20) for < x < 6 and 
a; < 0. 

4.4 Exact solutions for 0i in two special cases 

An exact analytic expression - derived from the general method in Section 



4.3 - exists for the screened potential (^1 in two cases. The first one cor- 
responds to Aci = so that the wall has no electrostatic response for any 
strength e of the Coulomb coupling inside the fluid. The second situation 
is the strict limit e = for Ad 7^ 0. Indeed, in both situations Ui = and 
when x' > b (pi is the solution (j)^^\x, x' , q; b, A^i) of 



dx'^ ^ 



dx"^ 



(/)(°)(x,i',q;6,Aei) = -47r5(x-f') if 5 < x (4.27a) 
0(°)(x,x ,q;5,Aei) = ifx<6 (4.27b) 



In the following we use the notation b for ki6 and introduce 5j only when 
we consider both Kib and K26. 

When Aei = 0, for any coupling constant e 

^1 (x, x', q; e, 6, Aei = 0) = (^(°) (x, x', q; 6, Aei = 0) (4.28) 



and has no special property at x = 0. (According to ([4. 5]) , (j)i and d(j)i/dx 
are continuous at the crossing of the wall.) In fact, up to a translation of 
its argument x equal to 6, 0i is defined as the mean-field potential in 
the linearized Poisson-Boltzmann approximation in the case of a multicom- 
ponent plasma in the vicinity of a plain hard wall (Agi = 0) located at 
X = 0. Such a wall exerts only a geometric constraint without any electro- 
static attraction : there is no need for introducing the hard-core repulsion 



Vsn (2^) involved in our generic model, so that b could be set to zero in 
this particular case. (We recall that at leading order in e, the correlation 
(or Ursell) function of charges with a hard-core diameter a is just equal to 
— /3eQ,eQ,'i;^(x, x'; y — y'; A^ = 0) |Q] and does not involve a.) As recalled in 
( 4.261 ) (pi may be written as the sum 



(/>(°)(x,x',Kiq;6, Aei = 0) =(I)^{ki\x - x'|,q) 

+ h%w{^ + x' - 26; Kiq) 
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(see Ref. ^]) where 0b is a particular solution of (4.27) which is chosen 
to be its bulk value, 



(/)b(|x - x'|,q) 



27r 



(4.30) 



and /i^tv is a solution of the associated homogeneous equation which van- 
ishes when X goes to +oo and whose coefficient is entirely determined by 
the boundary conditions at the interface, 



htwi^ + x -26; q) 



-{x+x'-2b)^/T+q^ 



(4.31) 



However (/'2(Aei = 0) cannot be calculated explicitly, because U2{x) ^ 
according to ( |4.13D and ( |4.29| ). 

When Aei 7^ 0, the electrostatic response of the wall disappears in the 
limit e = for 6 fixed, 



lim(/)i(x,x',q;e, 6, Aei) 

£— >0 



(x,x',q;6, A 



elj 



(4.32) 



where 



(?:)(°)(i,x',q;5, Aci) = (/>b(|S - i'|,q) 

+ Z(q;6,Aci)/i+^(x + x'-26;q) (4.33) 

Z is & renormalization factor arising from the continuity conditions at x = b 
when Aei 7^ 0, 



Z(q;6,A 



clj 



i-Aeie-25'' vr+7+iqi 



1 - A,,e-2gfc 



(4.34) 



4.5 General equivalent integral equations 

According to section the determination of a screened potential is equiv- 
alent to solving the homogeneous equation ( |4.21 ) . The asymptotic equation 
at large distances associated with ( [4.21| ) has two exact linearly indepen- 
dent solutions A± exp =pxyl + g^ which either vanishes or diverges when 

X tends to +00. Thus, in the following h^* will be looked for under the 
form 



hf*{x;q,b) 



(4.35) 



h^* is defined up to a multiplicative constant. For convenience sake, we 
choose the particular solution such that h^*(h; q, b) = exp (^b\/l + q^j ■ 
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When Uj{x) vanishes at large distances at least as fast as we show 
in Appendix ^that H^*{x) is the only solution of the integral equation 



where Cu- is a linear operator operating on a function / as 



(x) 



dv e 



■yv 



(4.36) 



(4.37) 



As it is the case for Dyson equation, the solution H^* of ( 4.36| ), denoted by 
Hij.^ in the following, can be written as the formal series 

H^* = H,.^^-C,.^ [1] + [1]] - [1]] ] + • • • (4.38) 

For the sake of conciseness we have omitted the dependence upon the pa- 
rameters 7 = 2^Jl + g2 and h in ( [4.38| ). 

We notice that, similarly, the function /ij"* which explodes exponentially 

fast when x tends to +00 and which is equal to exp[6y^l + q^] at x = h 
corresponds to a H^*{u;q) which obeys the integral equation 



dt e 



X l + Hr*{t;y/l + q^,b)\ (4.39) 
A series similar to ( 4.3^ ) can be written for H^* . 

4.6 Structure of the e-expansions of screened potentials 

The structure of the e-expansions of (pi and 02 can be investigated by means 
of the series representation (4.38) for the intermediate object Hu. combined 
with bounds upon Uy First we consider Ui{x) defined in ( 4.12|) , 



Ui{x) = y^ui< 



exp 



Acl£c 

2x 



1 



(4.40) 



with Mia = 47r/?e^2;Q,/Kj. As shown in subsection |B.l 

|Aei 



\U^(X)\ < 



E 



UiaSa exp 



^ ^ Ui a^a 

a 

Acie 



- if Ael < 

X 



2b 



1 



if Ael > 



(4.41) 
(4.42) 



The result of the e-expansion for the solutions * of the homogeneous 
equation ( 4.21 ) proves to be very simple at first order in e, as shown in 
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Appendix B. The term of order e in the e-expansion of h^* coincides with 
the corresponding term in the e-expansion of the solution of the equation 
( 4.7a| ) where Ui is replaced by its linearized value C/|™ 



(4.43) 



(We stress that such a coincidence is no longer valid at higher orders in e.) 
More precisely : 



hf*{x;q,b) 



1 + 



dv e 



dt e 



+ Oexp(e') (4.44) 



where e is a symbolic notation for the dependence on the e^'s. The double 
integral involving {7}™(t;e) in (4.44) is proportional to e, while C'exp(e^) is 
equal to - possibly multiplied by some power of Ine - times a function 
f{x; e, b) which decays exponentially fast over the length-scale yl + g^ ^nd 
which is bounded for all rr > 6 by a function of b and e/b. The function 
f{x;e,b) may contain powers x^^^-* (with lim^^Qp{e) = 0) times an expo- 
nential of — xyT+g^, as it is the case for the exact solution of the equation 
(14211 ) where C/i(i) is replaced by Uf'^iS:) (see (WM))- 

We point out that the correction of order e in h^*{x; q, b) is a function 
of X only, whereas the higher-order corrections are functions which vary on 
both scales 1 and e : 



h+*{i) = /if ) + (x) + eh'i^^ix) + e2/i(2) - 



X 



X,- ] +0{e-') 



(4.45) 



This property ensures that the zeroth order term in the e-expansion of 
the derivative dhf* /dx depends only on the corresponding term in the e- 
expansion of hf*, whereas the first-order term involves the derivatives of 
both the first and second corrections to /i^*, 



dh+* _ dhf^{i) 



dx 



dx 



+ e 



dh^^^^ii) dh}{'^{x,u) 



(2)- 



dx 



+ 



du 



U=x/£ 



+ 0{£^) (4.46) 



According to ( [4.39| ), results similar to ( 4.44] ) and ( [4.45 ) also hold for hi*{x). 

Finally, we turn to the e-expansion for the screened potential 0i. 
is calculated from the /i^*'s by using the expression (4.26) valid for x > b 
together with the continuity relations (4^) involving the solutions (|4.2C| ) in 



the region x < b. As a consequence of (4.46|), at leading order in e, the 



continuity relation ( [4. 5] ) for the derivative involves only the /if^^'s 
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and their derivatives dh^'^^/dx . Therefore, the leading term c^^^^ in the 
e-expansion of (pi is entirely determined by the functions hf''^ ^ and hf''^ 
(and not by functions which appear at higher orders in the e-expansions of 
the /i*^'s) and it coincides with (p^^'^ given in ( [4.33| ), 

(/)i(r,r') = Ki(^(0)(Kir,Kir';Ki6,Aei) + KiC'cxp(e) (4.47) 
where Oexp (&) has the same meaning as in ( [4 .441 ) . As a consequence of ( [1.47D 

-l(3el [(Pi - M (r, r) = e„F(i; e, h) (4.48) 

where F{x;e,b) is a continuous function of x in the interval b < x < oo 
which decays exponentially fast (with possible multiplicative powers x^^^^) 
over a typical length 1 when x goes to +oo. F{x;e,b) is integrable when x 
goes to +00 even when e vanishes, and may be bounded as follows. 



F{x;e,b) < MF{e/b)g{x) Vx > 6 



(4.49) 



6, +CXD 



where g{x) is continuous and integrable in the interval 

^From the previous result, we get the structure of U2{x) given in ( [4.13| ). 
For X = K2X 



Uo(x) 



^eaF(x(Ki/K2);e,b) 



(4.50) 



with U2a = ^T^P&tiZa exp(eQ,)/K2. We show in Appendix ^ that, as it is the 
case for Ui, the first terms in the e-expansion of h^* are given by ( |4.44 ) 
where ?7}™(5;;e) is replaced by 



U2''{x;£,b) = y^U2aea F{xki/k2;£, b) 



(4.51) 



The double integral involving U2^{t; e, b) in ( 4.44 ) is not merely proportional 
to £. According to definitions ( |4.8| ) and ( [4.S| ) 

K2 = [1 + 0{e)] (4.52) 

and the properties of the double integral at stake imply that the first two 
terms in the e-expansion of K^* are given by (4.44) where is replaced 

by 



^^U2a£a F{x;e = 0,6) 



(4.53) 



Subsequently, all properties derived from (f4.44|) also hold for /ig . 
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Eventually, the main result of the previous perturbative analysis is that, 
at leading order in e, 0i and 02 coincide with the same function <f>^^^ apart 
from a scaling dependence upon either ki or K2 



1 ~ 1 

(j)j{x,x',Kjq;e,Kjb,Aci) = —(p^'^^Kjx, njx' ,q; njb, A^i) -\ Oexple) 



(4.54) 



In ( [4.54 ) Oexp(e) denotes a function of the variable x which satisfies two 



properties. First it decays exponentially fast over the scale kj in the sense 
that it falls off as exp(— Kj-xy^l + g^) times a function which may increase as 
xP(^) with lime^oP(£) = 0. Second, OexpCe) remains bounded by e - possibly 
multiplied by some power of In e - times a function of e/ (nib) oc /3e^/6 for all 
X > b, even if Oexp{£) may drastically vary over the scale enj ~ /3e^ <C kj. 

As shown in next section, because of the neutrality constraint ( p. 131) , the 
explicit value of the correction of order e in (j)2 happens not to appear in 
the first correction to the density profile. (Only some properties of it must 
be known in order to settle that it really does not contribute to the density 
profile at order e). The latter result may be viewed as a consequence of the 
following property, already used in Ref.[Q|. According to the first equation 
of the BGY hierarchy, the first-order correction induced by Coulomb inter- 
actions in the density profile is determined only by the potential drop 
and the Ursell function at leading order; besides the latter one is nothing 
but the screened potential 02 in the limit e = 0. We recall that lim^^o 02 
is drastically different from the bare long-ranged Coulomb potential v^, 
which would be a too crude approximation for the Ursell function : though 
lim£_^o 02 has the same amplitude /5e^ as v^, it is a short-ranged function 
with a characteristic scale k oc e/(/3e^). 

5 Scaling analysis in the weak-coupling limit 

In the weak-coupling regime of interest, the particle density Paix) at leading 
order in e is equal to its value p^^i^) ideal gas submitted to the external 
potential corresponding to the self-energy —Aeie'^/4:X of a charge in the 
presence of a wall with an electrostatic response, 

p^{x) = exp (Aei/5e2 /4x) (5.1) 

According to ( |1.8| ) and ([4.§|), the small dimensionless coupling parameter 
inside the Coulomb fluid may be chosen as 

Ea = ^Ki^el < 1 (5.2) 
while the coupling constant with the dielectric wall 

Mi (5,3) 
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can take any finite given value. 



5.1 Screened fugacities 



In the weak-coupling regime, a simple scaling analysis may be performed in 
the bulk as in the inhomogeneous situation near the wall. As a consequence 
of ( 4.54| ) {l/2)Pe'^{4>i — VB){r,r) scales as £a at leading order and takes the 
generic form 

-^/3ea {(t>i - Vb) (r; r) L{ki{x - b); Kib, Ad)} + Oexp (e^) (5.4) 

where the term of leading order is entirely determined by 
, Kir; Kib, Aei), 



L(x -b:b,A 



el) 



(2vr)2 



(/)(°^ (x, X, q; b, Agi) - (Ab(5, x, q) 



(5.5) 



L{x — b; b, Agi) has been studied in Section 3.2 of Paper I. It can be decom- 
posed into 



L{x - 6; 6, Aei) 



-A, 



2x 



+ L{x;b, A. 



cl) 



(5.6) 



where L{x — b; b, Ad) remains finite even when x = b = 0. Thus, according 

[1] 



to ( ]3.5| ), ( |5.4| ) and (p^), Za proves to read 



(x) = Za e^'" 0{x-b) exp 



A 



el- 



Ax 



-2k\x 



{ 1 - e„Z(Kix; Kib, Ad) + Oexp (e^) } (5.7) 



The expression (|3.7] ) of z^' differs from the expression (^^) of only 
by the replacement of (pi by (f)2- At leading order in e, apart from the 
change of ki into K2, 02 coincides with 0i (see ( [4.54 )). Thus Za^ is given by 
(5.7) where ki is replaced by K2 and is multiplied by K2/K1. However, 
according to (|49|), 



= £«[! + 0{e)] 

Kl 



(5.8) 



so that 



(x — 6)zq, exp 



A, 



^ Ax 



{1 + Ea [1 - L{k2X; K2b, Ac 



+Oexp (e') } (5.9) 

We remind the reader that Oexp (£^) denotes a function which tends expo- 
nentially fast to a constant of order - possibly multiplied by some power 
of In e - when x goes to -|-oo. 
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5.2 Diagrams contributing at leading order in e 



First, we notice that even in the case where L is not bounded at the origin - 
as it is the case when ew > 1 ~ an e-expansion can be performed for integrals 
of the form j dv''Zly\x') f{r,r'). For that pmpose, we introduce 



wo{x;ea,Aei) = exp 



A ^" -2x 

Aei — e 
. 2x 



1 



(5.10) 



and we use the following formulas. Let us consider a dimensionless function 
/ of the variable x which involves two length scales h and l2- f may be 
written as fix/h^h/h)- We set e = h/h- When b > h, for all x > 5 x/h 
remains finite when e vanishes and 



Exp 

e^O 



dx f{x/l2;e) 



b/h 



dx2 Exp[/(x2;e)] 



(5.11) 



When b <^ I2, we use the identity 



Exp 

h^h Jb 



dx . . 



Exp 



[ dx ...+ [ 

Jb Jl 



00 

+ I dx 



(5.12) 



When I/I2 vanishes, X/I2 goes to zero for all x < I but remains finite for any 
X > I. Thus a basic formula for e-expansions when b <^ I2 reads 



Exp 

£^0 



dx f{x/l2\e) 



Exp Exp < li 

l/h^+oo 1/12-^0 [ 

+12 



l/h 

b/h 
+00 



dxi Exp [/(exi;e)] 



dx2 Exp[/(x2;e)] > 

l/h e^O j 

(5.13) 

We notice that only the sum of the two integrals in ( p. 13 ) is independent 
from /, whereas the first (second) integral may diverge when l/li becomes 
infinite (when I/I2 vanishes). Whatever the value of b may be, (5.11) and 
( 5.13| ) imply that Ku f dr' wq{x'] Eq, Aei) /(r, r') is at least of order elne (see 
also subsection 3.3 of Paper I). 

The simplest diagrams contributing to Pa{x) are written in (|3ll| ). The 
potential involved in F^^ is (f)2- First, the scaling property of (f)2 at leading 
order in e, namely 



4°^(r,r';K2) 



/t2</'*-°^(K2r, K2y') 



implies that, according to (|3J), (|4g|), (jSJ) and (^), 
j dr'z^^\x')F^''{r;r') = 0{e^) 



(5.14) 



(5.15) 
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Moreover, this leading correction depends on the species only through a 
coefficient z^. Subsequently, according to the neutrality constraint ( |2.13|) , 
after summation over species, the diagram with one bond F'^^ contributes 
only at next order in e. 



I dv'4^\x')F'"'{v;v')=0{e) 

7 

In the same way, the scaling argument shows that 
j dv'4^\x') [F-(r;r')]' = 0(e) 



(5.16) 



(5.17) 



An e-expansion can be performed for the contribution from [F'^'^] to the 
density representatio n (|3.11| ). Since L^kjx) decays exponentially fast, com- 
parison of (^]^) and ( ^.91) implies with ( ^.8]) that 



z\^'^ (x) = (x) X exp 



' 2x ^ 



+ O [z^e^) (5.18) 



and a straightforward calculation leads to the result 



dr' 



z^^\{x')-z^^\x')\ [F-(r;r')]' = 0(e2) 



(5.19) 



Theb asic formulcLS ( j5.11[ ) and. ( |5 . 1 3| ) with l\ — l3C(yC^ and I2 = 1^2^ allow 
one to show that 



j dv'z'^^\x')F^^{r-v') = 0{i' 



(5.20) 



Indeed, according to ( p. 10 ) and ( 5.14 ), the scaling change r = t/k2 shows 
that 



-^RT(f i r j £(17) 



If — f'l 



+ e„^</)2(f;f' 



K2a) e-^-^<^2(i^;f') - 1 
^:eLf^2(f;f'^ 



(5.21) 



with = K2/3eQ,e^. Since scales as Zq, = pj^^ oc 1/a^ (where p^^^ is the 
bulk density of an ideal gas with the same fugacities) and 



Exp FaT(f;f';ea7) = 0(47)' 



e-+0 



(5.22) 



/ dr'za'(x')-FRT(r;r') is of order ^2^0 ^^07 = ^i^lt-y)- (See definitions (|L 



and (|l.8D .) We notice that the same result is obtained by using the scal- 
ing change r = /?eQ,e^u and the property that lim^^o -^rt(£u, eu'; Gq,^) is 
an integrable function of u independent from e so that (|5.13D implies that 
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/dr'zL^'(x')FRT(r;rO is of order (/3e„e„/)3a-3 = 0(T^) = 0{£^). In fact, the 
lengths b and a may also generate contributions of order exp[Aei/3e^/46] 
or exp[/3e^/(T], where b and o" are supposed to be finite. 

A scaling analysis analogous to that performed in shows that more 
complicated diagrams contribute at least at the order in the exponential 
of the diagrammatic representation ( ^.11 ). Eventually, in order to calculate 
Pa{x) up to order e, in ( 3.11] ) we only have to consider the contribution from 
the diagram with one bond F'^'^ up to order e. 



5.3 £-expansions of the contribution from the F^'^ diagram 

In ( |3.11| ) the contribution from F''^^ up to order e is derived from ( |3.8| ), 
where (j)2 is replaced by its leading value 4>"2^ and (|5.9D . It is equal to the 
first order term in the e-expansion of 



r+co ^ 

J dx' z^e^ [l + wq{x; e-^, Aei)] 



X, x', q = 0; + (/.(I) (x, x', q = 0; 6, Aei) + axp(e')} (5.23) 



where b = K2b and wo{x; Eq, Agj) is defined in (|5.1[l| ). According to the neu- 
trality constraint ( 2. 131 ) and since the function (j) does not depend on the 
species 7, while the integral of /(x) times ^i;o(x; e^. Ad) gives a contribu- 
tion of order larger than / dxf{x), the contribution ( ^.23| ) involves only the 
leading term in the potential </)2; it is reduced to 



^2 Jh 



j dx' z^e^ |ti;o(x'; e-y, Aei) + 1 — L(x'; 6, Ad) | 



X 0(°)(x,5',q= 0;6) (5.24) 



where (/.W(i,x',q = 0; 6) = 27r [e"!*"*'! + e-(*+*'-2^)]. 
Equation ( 5.24| ) may be rewritten as 



j dr'z^^^{x')F''ir;r') = 

47r/3ea ^^g^^^ j]^ _ j^^^^-^^x; e^, Aei)]^ +C'exp(e^) (5.25) 



K2 



where = M + [M^ - M] with 



M(x;6,Aei) = ^ 



+ e" 



-{x+u'-2b) 



L(n';6,Aei) (5.26) 
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and — M is the e-expansion at orders Ine^ and (e-y)'^ of the integral 



1 1 

2 % Jb 



du' 



^~\x-u'\ ^ ^-(x+u'-2h) 



(5.27) 



5.4 Structure of the e-expansion for the density profile 

The structure of the e-expansion for the density profile can now be investi- 
gated. According to ( |3.11| ) and the results of the previous section 



where ^^'(x) is given by ( |5.9| ) and F'^"^(r, r'; a, 7) is replaced by — /5eae^K2 



1 + dr'z^^\x')F^^{r,v'-a,^) 



[l + Oexp(e')] 

(5.28) 



0(°)(K2r, K2r'). By using ( 5.25| ) we get 



p^{x) =e{x - b) exp ( 



4a; 



^ 1 1 + [1 - L (k2X; K2^, Aei)] 



z^e^e^ [1 - My (K2a;; e^, K26, Aei)] + C'cxp(e^ 



(5.29) 



By definition of the bulk density ( ^.11| ) and since L and vanish when 
X goes to infinity 



0(£2 



(5.30) 



The latter expression does coincide with the classical limit of the parti- 
cle density in a quantum plasma where exchange effects are neglected (see 



Eq.(5.28) in Ref.|Q]). Moreover, according to ( [Lq ) and (|4.9|) , K2 is equal to 
K]j up to a correction of order e 

= Kb [1 + 0(e)] (5.31) 
Using the property exp[— (1 + e)u] = exp[— — etiexp[— u] + O (e^u), com- 
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parison of ( p. 29 ) with ( |5.30| ) leads to 

=ple{x - b) exp (^Aei^e-^'^^^) 



X < 



1 - -^/3kd 



+Oexp(e^) 



(5.32) 



where = (l/2)/3e^K£,. 

The term with My in ( ^.32 ) is related to the electrostatic potential drop 
$(x) created by the charge density profile. Indeed, when the relation be- 
tween and the corresponding electrostatic field is combined with Gauss 
theorem together with the symmetries of the problem, and after an integra- 
tion by parts, proves to read 



r+co 

^{x) = -4tt dx' {x' - x)'^eaPa{x') 

J X 



(5.33) 



(The value of ^{x) when x goes to +oo is chosen to be equal to zero). When 
the structure ( ^.28|) of Pq,(x) is inserted into ( ^.331) , according to (|5.9|) (^.11 ) 
and ( |5.13| ), the contribution to (^{x) at leading order e only comes from the 
following part in Pa{x) (with x > 6), 



zPl(x) + z„(-/3e«)G(x) 



where 



G{x) = [dr'Y^ eyz^? (^') «2<^^°^ (/^2r, Ksr') 
Since (p^^'^ {x , x' , obeys equation ( [4.27a ), 



(5.34) 



(5.35) 



G{x) = 47r ^ e^z^^l (x) + 



d^Gjx) 
dx"^ 



(5.36) 



when X > b. Moreover k| = 47r/3 (J2a ^a^a) x [1 + 0(e)] so that 



$(x) 



+ 00 



dx'{x' — x) 



d^G{x') 



dx 



[l + Oexp(e)] 



(5.37) 
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According to the decay property of (j)^^'^ , G{x) tends to a constant exponen- 
tially fast when x goes to +00. Thus, after integration by parts, (|5.37| ) leads 
to 



(^ix) 



G{x) - lim G{x) 



[l + Oexp(e)] 



(5.38) 



Since —j3eaG{x) coincides with the l.h.s. of ( 5.25| ) at leading order in e, the 
explicit value of <I>(x) is given by 



(^ix) 



2^(3 

Kb 



6,Aei) + Oexp(£V/?e) (5.39) 



The latter result ensures that the expression of the density profile derived 
in Paper I from a mean- field approach is exact at leading order in e^,. 



A Appendix A 



In this appendix we show that, when U-^{x) vanishes at large distances at 
least as fast as l/i, H^*{x) defined in ( [4.35| ) is the only solution of the 

integral equation ( [4.36| ). H^*{x] \/\ + q^, b) obeys the nonhomogeneous 
differential equation 



and satisfies the following boundary conditions 



U^{x)H+ = U^{x) 



and 



lim e-^'^^H+^x; ^/l + q^,b) = 

x— >+oo 



(A.l) 



(A.2) 



(A.3) 



By using the extra change of function = exp —2xy'l + q^ 

xdH^ /dx the problem is reduced to a first-order differential equation which 
relates dGi/dx to Hj~. The formal integration of the latter equation leads 
to the relation 



H+{x;y^TT^,b) = ao + a^e^^V^ -Cu, \l + H+;2y^TT^,b\ (A.4) 



where Cu- is the linear operator defined in ( [4.37| ). 

When Uj tends to zero at least as fast as Uq/x when x goes to infinity 
(with Uq a constant), the boundary conditions ( [A .2] ) and ( |A.3| ) imply that 
the integration constants oq and ai are in fact equal to zero in the case of 
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H^* . Indeed, the hypothesis about U-^ means that there exists some xq 
such that 



Vx > xq \U\{x)\ < 



(A.5) 



On the other hand, according to ( |A.3| ), for any given real M > there exists 
some xi such that 



Vx > xi{M) e 



l + H^*{x)\< 



M 



(A.6) 



Since exp(7u) is an increasing function of v, ( A.5 ) and ( |A.6| ) combined with 
properties of integrals lead to the inequality 



l + H+*■2^/l + q^,b 



X2 



+ MUoe 



1; y^T+q^,X2 



-.2 



(x) (A.7) 



with X2 = sup{xo, xi). In ( [A.7 ), the left term in the upper bound is indepen- 
dent from X while the right term can be calculated explicitly. An integration 
by parts leads to the result 







1;7,6 


(i) = i " 


7 L 



ln(7x) + B{'jx; 76) 



(A.^ 



where B{^x;^b) is a continuous function of x which has a finite limit when 
X goes to +00 and which tends to — In 5 when x approaches the value b. 
Indeed, 



B{jx]jb) = - [e^*Ei(-7i) - e^^Ei(-76)l - ln(76) 



(A.9) 



where Ei(— n) denotes the Exponential-Integral function defined for u > 
as Ei(— ti) = — f^°° dt exp{—t)/t. As a consequence of ( [A.7 ) and (A. 8), 

in the large-x limit exp(— 2xy^1 + q'^)C,j. 1 + vanishes as well as 

X exp —2xy^l + q^ (according to ( [A.3[ )) so that ai in ( [A.4[ ) proves to be 

zero. Then the definition ( [4.37[ ) together with the choice ( [A.2| ) enforce that 
gq is also equal to zero. 



B Appendix B 

In the present appendix we study the structure of the e-expansions of the 
functions H^ defined by the series representation ( [4.38| ). 
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B.l Bounds for the i^rr 's 



According to integration properties, the linear operator Ly defined in (f4.37 ) 
has the following properties, 



l^c.[/]l<^|a|[l/l] 



(B.l) 



|C/|<C/'^£m[|/|]<£^4|/|] 



and 



0< /<5 



(B.2) 



(B.3) 



Therefore, by a recurrence argument, the series representation ( ^.381) of 
leads to the following result 



if \U\ < U' then \Hj,\ < H_jj, 



(B.4) 



In the present appendix we will omit all irrelevant indices and coefficients 
and we consider only the prototype functions 



Ui{x;e) 



1 



and 



C/2(x;e;6) = e"^(*'"'^)-l 



(B.5) 



(B.6) 



where F{x;£,b) is a continuous function of x in the interval b < x < +00, 
which is bounded for all x > 6 as follows. 



F{x;e,b) <MF{e/b)g{i) 

where g{x) is a continuous function integrable in the interval 
sequently for all ii > 6 

F{x]e~b) < MF{e/b)Mg{b) = M2{e/b~b) 



b, +00 



(B.7) 



. Sub- 



where Mg{b) denotes the upper bound of g in the range b < x < +00. 

Bounds upon the potentials f/j 's are derived by noticing that for all reals 

V 



le-l"! - 1| < \v\ 
< el"! - 1 < \v\e 



(B.9) 



( (|B.9| ) is a direct consequence of the series representation of the exponential.) 
(|B.9|) implies that 



(B.IO) 



29 



where ei = eMi{e/b, A^i) and 

Mi(e/6,Aei) = ^exp(^^0 if Aei > (B.lla) 
= ^ if Aei<0 (B.llb) 



Besides, if the sign of F{x) depends on x, F{x) may be written as F{x) = 
- with = F{x) if F{x) > and = otherwise, 

while -F-(x) = —F{x) if F{x) < and = otherwise. With these 

definitions = F^{x) + F^{x). Using the decomposition 

e'^ -1= {e'^+ - 1) e-'^~ + (e'^^" - l) (B.12) 
together with ( p^ ) and ( p^ , we get 

|t/2(i;e,6)| < f/2"'" = e25(i) (B.13) 
where £2 = eM2(e/6, 6) and 

M2(e/&, 6) = MF{e/b) exp [eAfF(e/^)M3(6)j (B.14) 
The bounds ( P-lOj ) and ( p.l3| ) about the C7j's ensure that, according to ( p.4| ) 

\Hu,\<H^,^H and \Hu,\<H^e,g (B.15) 

As shown in the following, some information about the structure of the e- 
expansions of the Hu.^s can be derived from ( B.15| ). In particular in both 



cases, the linearity of integrals ensures that the n^'^ term in the series ( [4.38 ) 



for i/_[/!nax is exactly of order e", which is not the case when U{x,e) is not 
a linear function of e. 

B.2 £-expansions of bounds for the Hjj.''s 

First we consider H-i^g. A mere integration by parts allows one to show that 
when git) is positive, integrable at large t and continuous for 6 < i < +00, 

/■+00 

Cg[l]{x)< y__ dtg{t) = G{b) (B.16) 

Since = — a recurrence argument proves that the n^'^ term in 

the series representation (|4.3j ) of -ff-eg is a positive function of x which is 

r _ 1 n 

lower than 

can be resummed with the result 

eG{h) 
1 - eG(b) 



£G(b) . Thus, H-sg{x) is lower than a geometric series which 
result 

H-sg{x) < -r^^ (B.17) 
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Therefore H^i;g{x; e,b) can be written as an e-expansion whose coefficients 
are bounded functions of x, 



+00 

E 

n=l 



(B.18) 



with Gn{x; h) < Gn{b) for all x >b. 

When U{x) = e/x the starting equation ( |A.1| ) satisfied by Hu, is exactly 
solvable and the formal series (|4.38 ) must coincide with the e-expansion of 
the explicit exact solution. More precisely, exp(— x) [l + H^/f{x)] is a solu- 
tion of the stationary radial Schroedinger equation for a quantum state with 
zero angular kinetic momentum and a negative energy -1, in the Coulomb 
potential —e/r in dimensionless units, and with boundary conditions ( |A.2| ) 
and ( [A.3| ). ( ( |A.2| ) is different from the corresponding condition in a Hydro- 
gen atom.) The equation reads 



dx"^ 



+ 



£-1 



/i = 



(B.19) 



The solutions of ( |B.19| ) are well-known : with the boundary conditions ( IaT 
and (lAlD 



1 + H,n{x) = ^(£,6)e^-%,/2,i/2(2x) (B.20) 



where A{£, h) 



a normalization factor 

W^./2,i/2(2x)cx(2xr/V^ 



1/2(2^) and W is the Whittaker function. Up to 

s/2 



1 + 



t 

2i 



(B.21) 



The structure of the e-expansion of H^ji can be derived either from the 
e-expansion of the above integral representation of We/2 , 1/2(22;) or from the 
series representation ( 4.3S ). ( ([4.38| ) is also an e-expansion of H^ji, since 
the n*^ term is exactly of order e" because of the linearity property of 
integration.) We choose to exhibit the second derivation in order to il- 
lustrate how the series representation provides information. A recurrence 
scheme using an integration by parts and the fact that £1/^/(7^); 7, h]{x) = 
(l/7)i2i/t[/(t); l,76](7x) allows one to show a generalization of ( |A.8|) (with 
7 = 2), 



i/i 



(ln2t)^';2,6 (x) 



^^-^ (In 2x)P+^ + (In 2x)" + i?,(x, 6) 

(B.22) 
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where Bp{x,b) is a bounded function of x when x > b. Since the n**^ term 



in the series ( 4.38 ) for U = e/t is exactly proportional to e", a recurrence 
scheme readily shows that the structure of the e-expansion of -^^/^(x) reads 



e/t( 



n=l 



1 /ln2x 



n! 



+ bn-i (In 2if-' + ... + bi ln(2x) + ^^(x, b) 

(B.23) 



In ( [B.23| ) B^{x,b) is a function which is bounded for x >b. The coefficient 
at every order vanishes when x = b (because £i/J/](x = 6) = 0). It 
diverges as (1/n!) (ln(2x)/2)" at large x and the summation over the leading 
asymptotic behaviours at all orders in e can be performed explicitly with 
the result 



~ e2 

->+oo 



ln{2x) 



(2x) 



(B.24) 



The asymptotic behaviour ( p.24| ) indeed coincides with the leading large- 
distance behaviour of the explicit solution (B.20) derived by using the inte- 
gral representation (B.21), 



e/2,l 



/2(25) 



(2x 



~\s/2 



8i 



(B.25) 



B.3 Exact solution at first order in e 



First, the fact that the upper bounds U^^^{x;£j,bys upon the C/j's are of 
the form 



?7j'"^"(x;ej,6) =ejl^'^^"(S) 



(B.26) 



(where ei and €2 are defined in (|b1o|) and (|b13|) ) allows one to show that 

+00 

Hu, = -Cu, [1] + Yl e"«j,n(*; b) (B.27) 

n=2 

with for all x >b, 

ai,niS:;e,b) < Mj{e/b,b) A^^nix,!^'^) (B.28) 



where the Mj's are defined in QB.llD and (|B.14D . In ( |B.28|) Ai_„(x; e, 6) does 
not diverge faster than (Inx)" when x goes to infinity, while A2^n{x; b) is a 
bounded function of x for all x > b. The demonstration is the following. 
According to the series representation ( [4.38| ) for Hu. , the inequalities ( [B.l| ) 
and (|B.2| ) together with the proportionality relation ( |B.26| ) lead to 



n=2 



(B.29) 
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where the coefficient of contains n operators iZ^max. ( B.29D imphes that 



the sum in ( B.27 ) indeed starts at the order e^. Moreover, according to 
( 4.38| ), the right term in ( |B.29| ) may be written as 



_ff_j^max -)- i2_ymax 
j j 



[1] 



(B.30) 



The exphcit structm'es of the e^-expansions (|B.18| ) and (B.23) for the 
ii'_(^max's together with the results ( A.8 ) and (|B.16D for £_[/max imply that 



+00 



(B.31) 



n=2 



where the have the properties given after (^]2|). 

Second, we show that the e-expansion of Lu. [1] — £rriin[l] starts at order 

^ j 

e^, where the linearized potential C/j^'" is the first-order term in the expansion 
of the exponential involved in the definitions (B.5) and ( |B.6| ) of the C/j's. 
More precisely 



£^.[1] =£^,i„[l] +e2i?^.(5;;e,6) 



(B.32) 



where for all x > b, R^. is a bounded function of x and the dependence of 
its bound ML.{e/b,b) upon e is entirely contained in Mi{e/b) if j = 1 or in 
[MFie/b)fexp[eMFie/b)Mg{b)] if j = 2. is derived as follows. First 

we notice that 



^Uj[f] - 'Cj,lm[l] = £y lin[l] 
J 3 



Moreover for all reals v 



l + \v\ \ <v^ 
,1 < i,2el«| 



so that 



>l _ 1 



C/i(f;e) - Ul'-{x;e) < \MM,{e/b)e'^ 

2 



(B.33) 



(B.34) 



(B.35) 



A decomposition similar to ( [B.12| ) reads 



l-eF= {e'^+ - 1 - eF+) e"^^- + eF+ ( 

+ (e 



-eF_ 



-1) 

1 + eF_) (B.36) 



and we get 



U2{x;e) - U^'''{x;e) < MF{e/b)M2{e/b,b)e^g\S: 



(B.37) 
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According to ( |B.2| ), the latter bounds together with ( p.33| ) imply that 



/:..iin[i] 



9j 



[1] 



(B.38) 



where gi{x) = {\A^i\/2)Mi{e/b)l/x^ and g2{x) = MF{e/h)M2{e/b~b)g{if. 
Since as well as g{x)'^ are continuous for x >b and integrable when x 
goes to oo, the result (|B.16D can be applied to the gj's and leads to ( B.32|) . 
Subsequently, according to ( [B.27] ) and ( B.32| ) 



Hr 



. (x) = -£^,in [1] (x) + e^Ru^ {x; e, b) 



(B.39) 



where if e/b and b are kept fixed, lim^^Q Ru.{x,£,b) < +00 and Ru.{x,e,b) 
does not diverge faster than x'P^^^ with lim£^Qp{e) = when x goes to +cxd. 



We recall that the result (B.39) holds for a potential C/j(x;e) which may be 
not bounded for all x > 0. The only hypotheses about are that 

Uj{x;£) tends to zero at least as fast as e/x when x becomes larger than 1. 
The property (KM can be applied to C/j'™ and we get 



HuAx) - H„uJx) 



Ru- (x; e, b) — R^i\u (x; e, b) 



(B.40) 



The solutions of the homogeneous equations ( [A.!]) for C/j and f/j™ respec- 
tively coincide only at the first order in their e-expansions. 
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Figure 2: Example for rule after step-1 resummations. After ring sum- 
mations, and contribute to G*J^\ whereas the subdiagram §(7-"!) 
in G*j^^^ disappears and G*J^^^ contributes to another G* . As a result the 

value of G*J^^'^ is equal to an integral where all bonds are the same as in G*j^^^ 
but where all weights z have been replaced by z^^] minus the value of the in- 
tegral corresponding to G*J^^^ where the weights z have been replaced by z^^' 
only for points which are not the intermediate points of a ring subdiagram. 




Figure 3: Step-2 for resummations. gJ and G^^ are two G* '1 diagrams 

which correspond to the same prototype diagram P as the diagram 
in Fig. |2[ The brackets with an index T indicate that the diagrams G* '^1 
must be calculated with the weights and the subtraction (rule JV) displayed 
in Fig. I 
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-[2] 

z 




Figure 4: Diagram P corresponding to , and G^j after step-2 

resummation. The crenelated line is a bond Frt, while the double wavyline 
corresponds to a bond (1/2) [F^'^]^. The rule DI2 is illustrated in the present 
example. 




-[2] _ [1] 
Z - Z 




P 

Figure 5: Example for rule Diagram G* contributes to the value of 
diagram G* '^1 after step-1 resummations. The value of G* is determined 
by rule "Ji. After step-2 resummations G* t"*^' contributes to diagram P, whose 
value is given by rule A single wavy line denotes a bond F'^'^. 
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